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<^ ■ Abstract 



By making use of the duality transformation, gauge field correlators of the Abelian Higgs 



fSj '■ Model are studied in the London limit. The obtained results are in a good agreement with the 



dual Meissner scenario of confinement and with the Stochastic Model of QCD vacuum. The 
nontrivial contribution to the quartic correlator arising due to accounting for the finiteness of the 
O ■ coupling constant is discussed. 

o : 

2P ■ 1. Introduction 

J . It is commonly argued that on the phenomenological level, quark confinement in QCD could 
be explained in terms of the dual Meissner effect [1] (for a recent review see [2]). That is why it 
looks reasonable to study various properties of the Abelian Higgs Model (AHM) [3-8] (see Ref. 
[9] for a review), where this effect displays itself, as a simple model of confinement. In this way, 
in Ref. [4] the duality transformation proposed in Ref. [3] has been applied to the lattice version 
^ of AHM in the London limit in order to reformulate its partition function in terms of the integral 
over world-sheets of the Abrikosov-Nielsen-Olesen strings. Then the same reformulation has been 
performed in the continuum limit in Refs. [7] and [8] for AHM with and without monopoles, 
respectively. 

In the present paper, we shall apply the duality transformation to the derivation of the string 
representations for the generating functionals of correlators of the gauge field strength tensors and 
of the Higgs currents in the London limit of AHM. By virtue of the obtained representations, we 
shall then get the bilocal correlator of the dual field strength tensors. It will be first obtained in 
Section 2 in a direct way, i.e. from the generating functional for the strength tensors correlators, 
and then in Section 3 from the correlator of two Higgs currents, by making use of the equations 
of motion. The asymptotic behaviours of the obtained correlator occur to be in a good agreement 
with the present lattice data concerning such a correlator in QCD [10]. The latter one plays an 
essential role in the Stochastic Vacuum Model (SVM) of QCD [11] (for a review see Refs. [2], [12], 
and [13]). This result supports the original conjecture of 't Hooft and Mandelstam concerning the 
dual Meissner mechanism of confinement. 
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In Section 4, we shall consider AHM in the vicinity of the London limit, i.e. take into account 
the effects bringing about by the fmiteness of the coupling constant. In this way, we shall see that 
these effects give a nontrivial contribution to the quartic correlator of the field strength tensors. 

The main results of the paper will be discussed in the Conclusion. 

In the Appendix, we perform a duality transformation of the generating functional for the field 
strengthes in the London limit of AHM. 



2. String Representation for the Generating Functional of the Cor- 
relators of the Gauge Field Strengthes 

We shall start with the following expression for the generating functional of the dual gauge 
field strength tensors in AHM 



Z[S aP ] = J |$|T?|$|TM^ sing W cg -exp 



I-/ 



4 ^ 



$1 



v 2 f + 



(1) 



flsing. + ^rcg n 



= 8^ — ieAp is the covariant 



where <&(#) = |$(x)| e t9 ^ is the Higgs field, 9 
derivative, and S a p is a source term. 

It is worth mentioning from the very beginning, that in what follows we shall be interested in 
the correlators of the dual field strengthes, rather than in the ordinary ones. This is because our 
main goal will be the derivation of the two coefficient functions D and Di, which parametrize the 
bilocal field strength correlator in the SVM (see Eqs. (12) and (13) below). Should we parametrize 
in such a way the correlator of the usual field strengthes, then as it has been explained in Refs. [2] 
and [12], the function D is nonvanishing only in the case of AHM with monopoles (otherwise, the 
term with the function D violates Abelian Bianchi identities). However, this problem is absent 
for the bilocal correlator of the dual field strengthes, in which the function D is nontrivial and by 
virtue of the equations of motion could be found from the correlator of two Higgs currents (see 
the next Section). 

In the London limit, A — > oo, one has |$| — > 77, and Eq. (1) takes the form 



[s aP ] = J vApVe^ve^expl-J d' 



x 



(2) 



Performing the duality transformation of Eq. (2) along the lines described in Refs. [3] and 
[5-7], we get 



[S af3 ] = JvApJ £>V J Vxp(i) exp I J d A x 



12r] 



- H 2 

2 n \iv\ 



mhp V Hp V - -F 2 „- 
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where E^x) = / da ill/ {x{^))8{x — x(£)) is the vorticity tensor current defined on the world-sheet 

E 

E of the closed Abrikosov-Nielsen-Olesen string, £ = ($, 1 ,^ 2 ) is a two-dimensional coordinate, and 
H^ u \ = + dxh^ + d v h\n is a field strength tensor of an antisymmetric tensor field h^ v . The 

details of derivation of Eq. (3) are presented in the Appendix. 
Now it is convenient to rewrite 



cxp 



d 4 xF, 



\1V 



/' 



VG fW exp\ I d 4 x -G^ + iF^G 



d x 



If 2 +iF 



2 hfw Sfiv 



after which ^-integration yields 

J VA^ exp 

= J VG^ exp (- J d 4 xG 

= J VB^ exp J d 4 x {^h^ + + d^B u - d u B 



y2 



(4) 



Next, performing the hypergauge transformation with the function -B^ (see e.g. Ref. [5]), we 
get from Eq. (4) 



Z [S a p] = exp ( - J d 4 xS j 
J 2*^(0 J Vh^exp | J d 4 x 



12rj 



1 e l 



(5) 



Gaussian integration over the field in Eq. (5) leads to the following final result for the 
string representation of the generating functional (1) 



Z [S a p] = exp (- J d 4 xS i 



■ J Vx^g) exp |- J d 4 x J d 4 y (7rE Ai ,(:r) + ieS x „{x)) D Xu ,p P (x - y) (irE^y) + ieS w (y))^ . (6) 

It is worth mentioning, that an analogous expression for the generating functional has been 
for the first time obtained in terms of another fields already in the last paper in Ref. [1] and 
further studied in the first three papers in Ref. [2]. There it has been done by making use of 
the techniques quite different from the ones applied in the present paper, i.e. without casting the 
generating functional into the form of the integral over the Kalb-Ramond field. However, the final 
results presented in these papers qualitatively correspond to our Eq. (6). 

In Eq. (6), the propagator of the field h^ v has the following form 

^A W W = < P (X)+< W (X), 

where 



3 
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^ + y (AT + K 2 ) J U A/ A P - J + 



2 b 



■ I 8 \pX u -\- 8 p V x\x p 8 p\x v x p 8 V pXpX\ 



(8) 



Here = Ki(m\x\),i = 0,1,2,3, stand for the modified Bessel functions, and m = erj is the 
mass of the gauge boson generated by the Higgs mechanism. It is worth noting that the term 

daxv(x) J dapp(y)D x 2 ^ p (x - y) 

E E 

could be rewritten as a boundary one, and therefore vanishes, since £ is closed. 

Let us now derive from the general form (6) of the generating functional for the correlators of 
the field strengthes the expression for the bilocal correlator. The result reads 



F x „(x)Fpp(y) 



1 



8 2 Z[S af3 ] 



Z[0] 8S Xl/ (x)8Spp(y) 



S„ (5=0 



= (8 X p8 up - 8 Xp 8p U ) 8{x-y) + 2e 2 



D x „,w(x -y)- 



-27r 2 /y da a p(z) J da y( (u)D a(3: pp(z -y)D^ M ( 



u — x) 



MO 



(9) 



where 



(...) 



-7T 2 / da a(3 (z) J da jC (u)D af3M (z - u) 

E E 



/2?a; M (0exp 



-7T 2 / d(T af3 (z) / rffT 7f (M)D a/3)7C (2; - U) 
E E 



and the term with the 5-function on the R.H.S. of Eq. (9) corresponds to the free contribution 
to the correlator. According to Eqs. (7) and (8), one can see that due to the large distance 
asymptotics of the modified Bessel functions, Dx Vt p P {x) has the order of magnitude e 2 r] 4 . Therefore, 
the contribution of the second term in the square brackets on the R.H.S. of Eq. (9) to the bilocal 
correlator of the dual gauge field strengthes is much smaller than the contribution arising from 
the first term, when 



erj 2 m < 1, 



(10) 



where |E| stands for the area of the surface S. Such an inequality holds at least in the weak 
coupling regime (which is in the line with the London limit) or in the case of sufficiently small rj. 
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Following the SVM, let us parametrize the bilocal correlator of the dual field strengthes by the 
two Lorentz structures as follows 



\D(x 2 ) + 



+- 



d 



A | Xp5 U p X p5 U p 



dv I x p S Xfl - x,,,8 



D 1 (V 



(11) 



In Eq. (11), the division into such Lorentz structures is chosen in such a way that one can due 
to the Stokes theorem rewrite the double surface integral of the structure with the function D\ 
as a boundary one, whereas this is impossible to do for the structure with the function D. Then 
in the approximation (10), by virtue of Eqs. (7) and (8), we arrive at the following values of the 
functions D and D\ 



m 3 Ki 
4n 2 \x\ 



and 



D 1 (V) 



m 



2tt 2 x 2 



k i rn / ' 
\xl + A K0 + K \ 



(12) 



(13) 



In Eq. (12), we have neglected the free 5-function type contribution. The asymptotic behaviours 
of the coefficient functions (12) and (13) at \x\ <^ — and \x\ ^> — read 



D 



D 1 



m 



Aty 2 x 2 



1 I |4 

ir z \x\ 



(14) 
(15) 



and 



D 



D 1 



m 



-m\x\ 



4V2tt 



m 



m \x\) 



3 ' 
2 



(16) 



(17) 



2\/27r2 ( m | x |)2 ' 
respectively. 

One can now see that according to the lattice data [10] the asymptotic behaviours (14)-(17) 
are very similar to the ones of the functions D and Di, which parametrize the gauge-invariant 
bilocal correlator of gluonic field strengthes in the SVM approach to QCD. In particular, at large 
distances both functions decay exponentially, and the function D is much larger then the function 
Di, whereas at small distances the function Di behaves as 7^ and is much larger than the function 

7 \x\ 

D in this limit, which is in the line with the SVM. This similarity in the large- and short distance 
asymptotic behaviours of the functions D and D\, which parametrize the bilocal correlator of the 
dual field strengthes in AHM and the gauge-invariant correlator in QCD, supports the original 
conjecture by 't Hooft and Mandelstam concerning the dual Meissner nature of confinement. 
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3. String Representation for the Generating Functional of the Higgs 
Currents Correlators 



In this Section, we shall rederive the coefficient function D in the bilocal correlator of the dual 
field strength tensors from the string representation for the generating functional of correlators of 
the Higgs currents. In the London limit, such a generating functional reads 



VA^VO^ Ve 1 ^- exp \ / d 4 x 



l -F 2 

4 » v 



rj- 



(d,9 - eA,) 2 + J^j, 



where = —ei] 2 (d^6 — eA^) is the Higgs current in this limit. 

Performing the same duality transformation as the one used in Section 2, we get the following 
string representation for Z [ 



Z [J^ = exp 



2 

m 


f d 4 xJ 2 (x) 


_~2~ j 





Vx^g) exp 



■exp|ee Al/Q/3 J d 4 xd 4 y -'- 







d 2 



dx a dy. 



tx j da a/3 (z) / da y( (u)D af3ilC (z - u) 



(x - y)\ Jp{x)J s {y)- 



+7rE jU/0 (y) \^-D x „,w(x - y) j Jp(x) 

Varying Eq. (18) twice w.r.t. J^, setting then J M equal to zero, and dividing the result by Z[0], 
we arrive at the following expression for the correlator of two Higgs currents 



(jp(x)ja(y)} = m 2 5 /3a 5(x -y) + e 2 e XvaP e^ pl(J 

d 



1 2 



+7T Z 



dcr S ((z) j da Xip {u)\—^D XvM (x - z) 



2 dx a dy^ 




D 



[x-y) + 



dy. 



-D 



{y - u) 



(19) 



It is easy to see that the contribution of the term (8) to the R.H.S. of Eq. (19) vanishes, 
whereas the contribution of Eq. (7) to the second term in the square brackets on the R.H.S. of 
Eq. (19) could be disregarded w.r.t. its contributon to the first term, if the inequality (10) holds. 
Within this approximation, making use of the equation^ 



{jp{x)j a {y)) 



d 2 



d 2 



\D((x-yf 



(20) 



y dx^dy^ dxpdy a t 

which follows from the equation of motion, we get from Eqs. (19), (20), and (7) the following 
value of the function D 



D(x 2 ) 



m 3 K\ (m \x\) 



4tt 2 



x 



(21) 



1 An analogous equation for the dual fields has been first presented in Ref. [2]. For the first time, the fact that 
the function D becomes nonvanishing in QED with magnetic monopoles has been mentioned in Ref. [12]. 
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where we have not again presented the free contribution. One can see that Eq. (21) coincides 
with the value of the function D given by Eq. (12) obtained from the string representation for 
the partition function of correlators of the dual gauge field strengthes. Notice also, that only the 
function D could be obtained from the correlator (19) due to the independence of the latter of 
the function D\. 



4. 1/ A- Correct ions to the Quartic Correlator of the Dual Field 
Strength Tensors 

In this Section, we shall demonstrate that accounting for the finiteness of the coupling constant 
A in AHM leads to a nontrivial contribution to the quartic correlator of the dual field strength 
tensors. Let us notice from the very beginning that we shall not be interested in free contributions 
to the correlators and therefore in what follows consider only the correlators of the fields defined 
in different space-time points. In the London limit, within the approximation (10), the dominant 
contribution to the correlator 



(F^ lUl (x 1 )F^ U2 (x2)F^ U:i (x 3 )F f , iUi (x^ , x± ^ x 2 ^ x 3 ^ x 4 , 

comes from the four-fold variation of the following term in the expansion of the generating func- 
tional (6) 



e 4 f 

— J d 4 xd A yd A zd 4 uS^(x)D^ Xp (x - y)S Xp (y)S a p(z)D a p r/(; (z - u)S l( (u). (22) 

The 4! = 24 terms containing all possible combinations of indices and arguments, which one gets 
during this variation, are obvious and we shall not list them here for shortness. 

Let us now consider AHM in the vicinity of the London limit. To this end, we shall expand 
the radial part of the Higgs field, |<&(x)|, in Eq. (1) as |3>(x)| = rj + Tip(x), where r = j- — > 0, and 
ip(x) is an arbitrary quantum fluctuation. Correspondingly, in what follows we shall keep only the 
terms linear in r. Then neglecting the Jacobian emerging during the change of the integration 
variables, |3>(x)| — > ip(x), which will be eventually cancelled after division of the final expression 
for the field correlator by Z[0], we arrive at the new expression for the generating functional (1), 
which reads 



Z [S a( j] = / V^VA^V6 smg -V6 ICS - exp\- d A x 



(23) 



Following the same steps which led from Eq. (1) to Eq. (6), we get from Eq. (23) an additional 
weight factor in the integral over string world-sheets standing in Eq. (6). This weight factor 
emerges due to the -^-integration and reads 



Ki (M \x - y\) 



\x 



(24) 



In Eq. (24), M = 2y2A^ is the mass of the Higgs field, and Hffi- stands for the strength tensor of 
the saddle-point value of the field h pu following from Eq. (5). Taking into account that d^Y,^ = 0, 
since £ is a closed surface, we get the following result for this saddle-point value 
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hTfix) = £- 2 J ^ /'^Y' 1 ^ [i^{y) - es xv {y)\ + 



+-d p [d„S Xp (y)-d x S up (y)}\. (25) 



e 



Let us prove that the terms with the derivatives of S^ u on the R.H.S. of Eq. (25) yield zero. 
Due to the Hodge decomposition theorem (see e.g. [14]), S^ u could be always represented in the 
form Sfa, = d^Ku^dyK^+daLa^, where and L afll/ stand for some vector and an antisymmetric 
rank-3 tensor respectively. In the product S^F^, the contribution of the vector obviously 
vanishes due to the partial integration, and we are left with the representation of S^ v of the form 
Sfj, v = d a L aflu . Substituting it into Eq. (25) we see that the terms with the derivatives of on 
the R.H.S. of this equation vanish. 

Finally, upon substitution of the rest of Eq. (25) into Eq. (24), we arrive at the following 
additional term in the expansion of the generating functional, whose contribution to the quartic 
correlator in the approximation (10) is dominant 

J d 4 xd%d\d 4 vd A zd 4 w Kl ( ^ ~ yl) D 1 ((x - uf) D 1 ((x - v) 2 ) D, ((y - zf) ■ 



■Di ((y - w) 2 ) (x - u)^(y - z) a 



[x - v)^(y - w) a S u x{u)S uX {v)S[] 1 (z)S[]~ f ( w ) + 



+4(x - v) x (y - w) 1 S uX (u)S lxu (v)Sp 7 (z)S a ^(w)+ 
+A(x - v)^(y - w) J S„ x (u)S„ x (v)Sp 1 (z)S a p(w) 



(26) 



In Eq. (26), D x stands for the coefficient function (13) entering the bilocal correlator in the 
London limit. Even without explicit writing down 3 • 4! = 72 terms following from Eq. (26) 
after its four-fold variation, we see that the resulting quartic correlator differs from the one of the 
London limit, which could be obtained from Eq. (22). The most crucial difference is due to the 
presence of the Higgs boson exchange in Eq. (26), which is absent in Eq. (22). 

The other important outcome of Eq. (26) is that the leading correction to the quartic correlator, 
which arises due to the finiteness of A, could be completely described in terms of the coefficient 
function (13), if the inequality (10) holds. 

5. Conclusion 

In the present paper, we have studied dual gauge field strength correlators in AHM both 
in the London limit and beyond. To this end in Section 2, the generating functional for these 
correlators in the London limit has been cast into the form (6) of the integral over world-sheets of 
the Abrikosov-Nielsen-Olesen strings. In approximation (10), this generating functional yielded 
the values of the two coefficient functions, D and D ± , which parametrize the bilocal correlator 
and play a key role in SVM of QCD. These functions are given by Eqs. (12) and (13), and 
there asymptotic behaviours (14)-(17) at small and large distances are found to be in a good 
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agreement with the known lattice data [10] concerning the corresponding behaviours in QCD. 
This fact together with SVM supports the conjecture of 't Hooft and Mandelstam about the dual 
Meissner nature of confinement. In Section 3, by making use of the equations of motion, we 
have rederived the coefficient function D from the string representation for the correlator of two 
Higgs currents (19). In Section 4, we have studied AHM in the vicinity of the London limit and 
demonstrated that there appears a nontrivial contribution to the quartic correlator of the dual 
field strengthes. This contribution is due to accounting for the finiteness of the Higgs boson mass, 
which leads to intermediate exchanges by this boson in the quartic correlator. Besides that, we 
have demonstrated that according to Eq. (26), in approximation (10), the quartic correlator near 
the London limit could be described only in terms of the function D\, given by Eq. (13). 
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Appendix. Derivation of Eq. (3). 

In this Appendix, we shall present some details of derivation of Eq. (3). Firstly, one can 
linearize the term {8,9 — eA,) 2 on the R.H.S. of Eq. (2) and carry out the integral over # reg ' as 
follows 



j V6^- exp i~t J d 4 x (8,6 - e^) 2 j 



J VC,V6 rcg - exp I J d 4 x 



2rj- 



;pi + iC li (d li 6-eA ll ) 



= j VC,5 (8,C,) exp | J d 4 x -±- 2 Cl + iC ^ (d^ - eA,) 



(A.1) 



The constraint 8,C, = could be uniquely resolved by representing C, in the form C, = 
^£fMv\pd v h\p, where h\ p stands for an antisymmetric tensor field. Notice, that the number of 
degrees of freedom during such a replacement is conserved, since both of the fields C, and h, v 
have three independent components. 

Now, taking into account that s pu \pd\d p 6 snig \x) = 2nT lpu (x) [3], we get from Eq. (A.l) 



/ 



V0 sins ve ieg - exp | -y / d 4 x (d p 9 - eA^f 



Vh pu Vx p (() exp < J d 4 x 



ze 



(A2) 
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In derivation of Eq. (A. 2), we have replaced X># smg - by (since the surface £ parametrized 

by is just the surface, at which the field 9 is singular) and for simplicity have not taken into 
account the Jacobian arising during such a change of the integration variable |. 

Finally, adding the term —iS^F^ to the argument of the exponent standing on the R.H.S. of 
Eq. (A. 2), we arrive at Eq. (3). 
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